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Abstract 

This paper presents the operator form of the effective potential V 
<D \ governing the time evolution in 2 and 3 and n dimensional subspace 

of states. The general formula for the n dimensional case is considered 
the starting point for the calculation of the explicit formulae for 2 and 
^ | 3 dimensional degenerate and non-degenerate cases. We relate the 2 

!-h ' and 3 dimensional cases to some physical systems which are currently 

investigated. 
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1 Introduction 

The time evolution of physical systems in the Hilbert space is described by 
the Shrodinger equation: 

ij t \il)-t>=H\4,-t >, (1) 
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where % = c = 1. 

If we choose the initial conditions: 



hM = o>=K>>, (2) 

then the time evolution is described by a unitary operator U(i)\ip >= t > 
(|V >,\ip;t>E H, U(t) = e~ itH ). 

Vector \i[);t >G TC carries complete information about the physical system 
considered. In particular, the properties of the system which are described 
by vectors belonging to a closed subspace TC\\, of TL can be extracted from 
\if}',t > . In such a case it is sufficient to know the component \tp;t >iiG TLu 
of \ifj',t >■ The subspace TC\\ is defined by a projector P: H\\ = VTC, which 
simply means that \ip; t >\\= P\ifi; t >. 

Alternatively, the same result can be obtained by studying the time evo- 
lution not in the total space of states Ti but in a closed subspace Ti\\. In 
this way the total state space is split into two orthogonal subspaces H\\ and 
H± = Ti 7i\\, and the Shrodinger equation can be replaced by equations 
describing each of the subspaces respectively. The equation for H\\ has the 
following form jl] — [[|: 

- PHP^j Hj; t >||= | X ; t > -i | o °° K{t - r)|^ r >„ dr, (3) 

Q = I-P, (4) 
K(t) = Q(t)PHQe~ iQHQ QHP, (5) 
\ X ;t>=PHQe- lQHQ \ij>^ (6) 



where 



1 for t > 
for t < 



e(t) = 

Of course K(t) ^ only if [P, H] ^ 0. Condition (f3) can now be rewritten 



as 



IV; t = o >n= |v >n, IV; * = >±= |V >±, (7) 
where |V >±= Q|V > • 

If we now assume that at the initial moment no states from Hj_ are occupied, 
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>j_= 0, ( that is \x',t >= 0, \ip >= -0 >||) and define the evolution 
operator for the subspace 7i\\: 

\^;t >||= PKM P27(t)|V >= Pi/(t)P|V >||, (8) 

so 

U [{ (t)\iP >^ PU(t)P\iP > lb (9) 
we can transform (|3|) into 

^-Ptfpj U\\\if> >,|= -i j™ K{t-T)U\\{r)\il) >|| dr. (10) 



An equivalent differential form of (|T0|) has been found by Krolikowski and 
Rzewuski [0, §: 

^-^||(*)V||WI^>||=0, *>0, f/|,(0) = P, (11) 

where the denotes the effective Hamiltonian: 

H^t) = PHP + V\\{t). (12) 



For every effective Hamiltonian H» governing the time evolution in Tin = PH., 
which in general can depend on time t I0J - , the following identity holds 

1-0: 

where [[/||(£)] _1 , is defined as follows 

[f/ii(tr i [/ii(t) = f/||(t)[^i(t)]- i = p. (M) 

In the nontrivial case 

[P,H]^0, (15) 
from (|13|), using fll2|) and (§) we find 

p,i(t) = p/j^p^)]-^ (16) 

= PiPP+M,(i). (17) 
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and thus 

V\\{t) = PHQU{t)[U\\{t)]- l P. 

Assumption (p~5f ) means that transitions of states from Tin into Tl± and from 
H± into TC\\, i.e., the decay and regeneration processes, are allowed. Thus || 

- §, 

H ]{ (0) = PHP, V||(0) = 0, V|,(t -> 0) ~ -itPHQHP, (18) 

so, in general H\\(0) ^ > t = 0) § — § and Vf| (t ^ 0) ^ Vj^t ^ 0), 
if||(t 7^ 0) 7^ ^ift* 7^ 0)- According to the ideas of the standard scattering 

theory, it can be stated that operator Hu(t — > oo) = iifii(oo) = f /fn describes 
the bounded or quasistationary states of the subsystem considered (and in 
this sense it is similar to e.g. the LOY-effective Hamiltonian [§). 

From (|TU|) and (0),dg) it follows that the action of Vf|(t) on U\\(t) has 
the following form: 

POO 

Vf,(t)E/j|(t) = -i K{t-r)U\\{r)dr. (19) 

J 

The approximate form of V\\ can be obtained from (|10D and (|T9|) with the use 
of the retarded solution of: 

^-PHP ) jG(t) = P5(t), (20) 

where G(t) is the retarded Green operator: 

G = G(t) = -t&(t)e- itPHP P. (21) 

Then, using the iteration procedure for the equation fllPf) for [/n [[|, ^ ||. 
we get: 

oo 

= E/J(£) + YX-i) n L o L o L o ... o L oUtijt). (22) 



n=l 



n times 



Uu(t) is the solution of the following "free" equation H || 0: 

^ - PIfp) t/jf(t) = 0, t/fl(0)=P. (23) 
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o stands for the convolution 



f°g(t)= I f(t-r)g(r), 

and 



o 



L = L(t) = G o K(t). 



Equations fll9D and (g^) yield: 



V\\{t)U\\{t) = -%K o V\(t) + J2(-i) n K 0L0L0...0L0 U°(t). (24) 

n=l 

If ||L(£)|| < 1 then the series ( |24] ) is convergent. It is worth noticing that, 
unlike in the standard perturbation series, it is not necessary for the pertur- 
bation H l to be small in relation to H° (the full Hamiltonian H = H° + H 1 ) 
if ||L(£)|| < 1. This is considered one of the advantages of this approach over 
the standard ones as it can describe both week and strong interactions ||. 
If for every t > ||£(£)|| < 1 then to the lowest order of L(t), Vu(t) is 
expressed by ||: 

V\\(t) ~ V[f(t) = -i K{t - r)e i ^ PHP Pdr. (25) 

This formula was used to compute V\\ (t) for one-dimensional subspace TC\\ and 
to find the matrix elements of Vn(t) acting in a two-state subspace Tin in [^]. 
In some problems it is more useful and more convenient to use the operator 
form of V\\ (t) rather than the the matrix elements of V» (t) only. Searching for 
the global transformation properties of Vfi (t) under some operators expressing 
symetries of the system is as an example of such problems. 

Result ( pop will be the starting point for the following considerations 
concerning the explicit operator form of V\\(t) in n, 2 and 3 dimensional 
cases. 



2 Effective potential V\\ in n-dimensional 
subspace of states. 

Let us consider a general case of effective potential Vj| jn , acting in an n- 
dimensional subspace of states. Formally, the equation corresponding to 
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Eq.(|T2"|) has the following form: 

H ln (t) = PHP + V ln (t). (26) 
The projector P is defined in the following way: 

p = Y,\*i ><^\ = h, ( 27 ) 

where I\\ is the unit operator in TC\\, {|e_j > }jeA an d {|ej > }j=i,2,...n C 
i\ej > }jeA are complete sets of orthonormal vectors < e^e* >= 5jk in TC 
and Tin C respectively. Consequently, if the state space for the problem is 
H then H\\ = VH and P is the unity in H\\, P = I\\ 0. 
The subspace TC\\ can also be spanned by the eigenvectors of the hermitian 
matrix PHP: 

PHP\\j >= \j\Xj >, w=i,2,..., n ). (28) 

Using \Xj > we define projectors Pj |J where for simplicity the non- 
degenerate case of Xj is assumed. : 

Pj = - Al - lAiXA.i (M 3). (29) 
Of course these projectors fulfill the following completness condition: 

n 

£ Pi = p (so) 

The operator PHP can now be written as follows: 

n 

PHP = Y, X iPji ( 31 ) 

and following: 

Pe ±itPHP = P ^ e ±^p. (32) 
i=i 

This result can be directly applied to equation (|25| ) 

Vf|,„(t) = -i£ /"* PHQe-^ ~ T )(Q H Q '- X j)QHP dr P r (33) 
i=i 
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The integration can be easily performed, with the result: 

n -it(QHQ-\j) _ i n 

V\Ut) = -iY,{PHQ nMn r QHP}p j — ^2E(Xj,t)Pj, 



where 



-it(QHQ-X) _ i 

B(A, «)gPgq QHQ - A ^ 



Knowing that |TD 



1 



lim S(A.t) = PHQ — 
*-oo v ' ; QHQ -A -HO 



and defining 



we finally get: 



^ = PHQ QHQ-X + l QHP 



(34) 



(35) 



(36) 



V ln ^ fim V hn {t) = -zhrn £5^,*)^ = -i^E(A i )P J , (37) 



i=i 



3 Vii in a two dimensional subspace. 



In this section we find the explicit formula for V\\ in a two-dimensional sub- 
space of states using the framework presented above. 

In this case PHP, being a [2 x 2] hermitian matrix, has the following form: 



PHP 



H\\ H 12 
H21 H22 



where Hj^ =< £j\H\ek >■ 

The eigenvalues of PHP are easy to calculate: 



PHP\\ >- 



H\\ H12 
H21 H22 



OL\ 

0(2 



Oil 
0(2 



A ' = -(H n + H22) ± W|#i 2 | 2 + ^(#11 - H22) 



(38) 



(39) 
(40) 
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and if we adopt the symbols used in 0-0: 



A 1 ' 2 d ^ f H ± k, 



where 



Ho = -^(Hu + H22), 



K 



\H 



12 



+ 4(^11 



22 1 



Following, the eigenvector |A X > can be chosen as follows: 

#12 \ 



IA 1 >= 



#o+k-#ii 
1 



and the projector Pi. 



Pi 



1 



< Ai|Ai > 

/ #12 



|Ai >< Ai 



|#1 2 | 2 



+ 1 



Ho+k-Hu 
1 



H 



21 



H Q + K- H n 



(H +K-H n )-\H; 
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|ifl 2 | 2 + (i?0+K-^ll) 

{H +k-H 11 )-H 21 



(ff +«-ffii)-ffi 2 



so, explicitly 

Pi = 

L \H 12 \ 2 +irr ll +H-ifnf \Hi2\ 2 ( ^,+/.- //,,)- 

For clarity let us define: 



-(fl-O+K-ffll) 5 

(go+ K -H n ) 2 



dcf 



Pll Pl2 
P21 P22 



= 1,2). 



Both Pi and PHP can be represented by Pauli matrices: 

Pi = Po CT o + pI^w + P l y (?y + pI(t s , 
PHP = H a + H x a x + H y a y + H z a z , 
and the calculation of the coefficients p> yields: 

PO = |(pn +P22) = i 

= |(pi2+P2l) = 2#, 
^ = |i(pi2-p 2 l) = 2^, 
^= |(P11-P22) = 



We can see from the above that p u , (y = 0, x, y, z) can be expressed by 
H u , [y = 0, x, y, z), so finally we get the following expression for Pi 

p ^l^-^h + l PHP \ ^ 

Keeping in mind the fact that in 7i\\ we have cro = In = P, we obtain: 

P ^\((^) P + \ PHP \ < 48 > 

and after performing the same calculation for P 2 . 

p ^W + f) p -\ pHp \ ^ 



It is easy to verify that the completness condition Q50D is fulfilled: 

Pi + P 2 = P 



If we now come back to Eq. (|34| ) and use the results obtained in this section, 
the effective potential Vfi will have the following form: 



V\\(t) = -~E(H + K,t) 



K K 



-~H(# - k, t) \(1 + **°)P - -PHP] . (50) 

Matrix elements of this V\\(t) are exactly the same as those obtained in [^]. 

As noted in J7[ this result is significant. For example in the case of neutral 
K mesons the assumption of CPT invariance and CP noninvariance in the 
quantun theory, that is [CPT, H] — and [CP, H] ^ 0, yields: 

/Hi - h 22 ^ 0, (51) 

where hij =< ei\H\\\ej > are the matrix elements of Hu = PHP + Vu, 

V\\ = f V[|(oo), which runs counter to the usual assumption. More remarks on 
this problem can be found in the conclusions. 

The case of both eigenvalues of PHP equal can easily be obtained from 
the general case described above. The assumption of both eigenvalues equal 
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for a hermitian [2 x 2] matrix yields Hn = H22 and Hyi = H21 = 0. It is 
easy to verify that A 1 = A 2 k, = 0. Then: 

A 1 = A 2 = H (52) 

PHP = H n P 



and 



p e itPHP = e itH 0p ( 53 ) 



Thus, from equations ( f25"|) and (|33|): 



V||(t) ^ T^i(t) = -S(J^ ,^)^ (54) 

Furthermore, if apart from assuming the degenerate case of PHP we take 
t — > 00 we will get the same result as obtained from the Wigner-Weisskopf 
approximation by e.g. Lee, Oehme and Yang ||. It is interesting to notice 
that in this case hn = /i 2 2 (where hjj =< j\H\\\j >) with [CPT,H] = 0, 
whereas in the case of A 1 7^ A 2 under the same conditions we have (|5l|). 



4 Mi in a three dimensional subspace. 

This section describes the explicit formula for Vn in a three dimensional 
subspace of states in a very similar way as it was done for the two dimensional 
case. 

In this case the PHP matrix is a [3 x 3] matrix, for example 





' H n 


H\2 


H13 


PHP = 


H21 


H22 


H23 




_ H 3 i 


H32 


H33 




Hij = 


TT* 

Hji, 





and has the following characteristic equation: 

A 3 + AX 2 + BX + C = 0, (56) 
A = -(H n +H 22 +H 3a ), 

B = H n H 2 2 + H n H 33 + H22H33 — \Hiz\ 2 — \H 2 3\ 2 — \Hi 3 \ 2 , 

C = —(H11H22H33 + 2Re(H 12 H23H 3 i) — Hn\H 2 3\ 2 — H 2 2\H n \ 2 — if 33 |if 12 | 2 ). 
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It is easy to notice that A, B, C G 3? so, given the fact that PHP is a 
hermitian matrix, equation (56) is a third order equation with real coefficients 
and real solutions. To find the solutions we will use the Cardano formulae. 
Bearing in mind that the solutions are real we get the following three cases 
Ai 7^ A 2 7^ A 3 , Ai = A 2 = A 7^ A 3 and Ai = A 2 = A 3 = A: Let us find the 
eigenvectors, projectors and the quasipotential for each of the above cases. 



4.1 A x ^ A 2 ^ A 3 . 

In this case the three solutions of the characteristic equation 
by the following formulae: 



are given 



Ai 
A 2 
A 3 



A 2 


-3B 




3 


A 2 


-3B 




3 


A 2 


-3B 


3 



1 2 cos |a 



A 
3 



) 2 cos I (a + 27r) 
12 cos 4(a + 4-71") 



3 • 
A 
3 ' 



(57) 



<-+C) 



where cosa= — 5—^ . 

3<- 3 1 2 

The following basis of orthogonal eigenvectors can be chosen: 



w >-- 



A,) 



H23H12 
H13H21 



\Hi2\ 2 — (H n — Xj)(H 2 2 — Aj) 



1,2,3. 



(58) 



where j 

Using these eigenvectors we create projectors P in the way given in Sec. 2.: 

1 



< A,|A„- > 



(59) 



{\Hi 3 (H 2 2 — Aj) — i/ 2 3-^12| 2 + 



\H 2 3(H 



11 



\H 



12 



x 



A, 
-(H n 



- H 13 H 2 i\ + 
XMHn - Xj)] 2 }- 1 



x 



PU Pl2 Pl3 
P21 P22 P23 

di P32 P33 



0=1,2,3), 
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where 



Pll = 


#13 (#22 


- A,) 


— H 23 H 12 2 , 






Pl2 = 


(#13 (#22 


-V 


) — #23#12)(#32(#11 


-a,; 


) — #31#12), 


pL = 


(#13 (#22 


-a,; 


) — #23#12)(|#12| 2 — 


(#11 


-A,)(# 22 -A,)) 


P21 = 


(H 23 (H n 


-a,: 


) — #l3#2l)(#3l(#22 


-a,: 


) — H 32 H 21 ), 


P22 = 


I #23 (#11 


- A,) 


— H Vi H 2 i 2 , 






P23 = 


(#23 (#11 


-a,; 


) — #i3#2l)(|#12| 2 — 


(#11 


- X 3 ){H 22 - X,)) 


P31 = 


(\Hu\ 2 - 


(#11 


— Xj)(H 22 — Xj)(H 31 


(#22 • 


- Xj) — H 32 H 21 ), 


P32 = 


(\Hi2\ 2 - 


(#11 


~ Aj)(#22 — Xj){H 32 


(#11 


- Xj) — H 3 iHi 2 ), 


P33 = 


(l#12| 2 - 


(#11 


- A i )(# 22 -A i )) 2 , 







(where j = l,2,3). 

And consequently the quasipotential 

V|| 1 3(t) = -*E s ( A i > *) p i. ( 6 °) 



4.2 Ai = A 2 = A ^ A 3 . 

In this case we have the following expressions for the solutions of the char- 
acteristic equation 



\ _ / 2A A B 1 r \\ A 

^ ~~27 JT ' "3 ' (Ql) 

A 3 = -2( 2 #-f + C)l-f (61) 

In this case to define one of the projectors, say P 3 we can use the result 

presented above, so the projector will be given by formula (|59|). 

We do not actually need to know the remaining two projectors explicitly as 

V l3 (t) = -iS(A, t){P x + P 2 ) + H(A 3j t)P 3 (62) 

and Pi + P 2 = P — P 3 , P is the unity in the considered space so: 

V l3 (t) = -i~(X,t)(P-P 3 ) + ~(X 3 ,t)P 3 (63) 
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4.3 Ai = A 2 = A 3 = A. 



This case is the simplest one, and the solutions are: 

X = -j = H n = #22 = #33 (64) 

In this case PHP is a diagonal matrix in any basis. In fact, this is true for 
any n-dimensional hermitian matrix with all eigenvalues equal, so we get a 
form of quasipotential which is identical to the two dimensional degenerate 
case ([54]) . 

Vf| = -E(A,t)P (65) 

Again, if apart from assuming the three-fold degenerate case of PHP we 
take t — > oo we will get a result which is analogous to the one obtained from 
the Wigner-Weisskopf approximation by e.g. Lee, Oehme and Yang M. 



5 Equation for p matrix in Tt\\ . 

This section contains one possible application of the result obtained above, 
which is the equation for the density matrix p in Tin. 



Very often systems of the type described in Section 1. are considered 
as open systems interacting with an unknown rest, i.e., with the reservoir 
11, 12]. Then, for the description of the time evolution in subspace Ti,u 



instead of the state vector \ip;t >\\ solving equations (|]), ([TT]) , density matrix 
p is used. The p-matrix in quantum mechanics fulfills the following equation: 

^P = i\p,lT\, (66) 

where H is the total Hamiltonian of the system under consideration acting 
in the Hilbert state space 7i. H and p are hermitian. 

The consideration of such systems sometimes begins with a phenomenological 
Hamiltonian H e ff = H», acting in an n- dimensional subspace TC». Such 
Hamiltonians are of the LOY type or the type used in the master equation 



approaches [11, 12]. These Hamiltonians are not hermitian, therefore the 



time evolution of the reduced p-matrix, i.e., p\\ (where p» denotes the part 
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of p-matrix acting in 7i\\ ), is given by [O] 



d 

ft"' 



hHwpw -p\\ H \\ 



(67) 



where 



/ Pn{t) Pn(t) 
Phit) P22(t) 



P\\(t) 



Pin® Plait) 




Plnit) 
P2n{t) 



Pnn{t) 





V 



(68) 



At this point one remark concerning the above should be made: all properties 
of p\\ (t > 0) solving this evolution equation are determined by the form and 
properties of H\\ so for the same initial conditions po but different H» a 
different p\\(t) can be obtained. 

Let us notice that the solution of Eq.(|5T)|) has the following form 



p(t) = U(t)p Q U + (t), 



(69) 



where po = p(0) and U(t) is the total unitary evolution operator for the 
system considered. From this we conclude that the exact reduced p-matrix 
for a given complete and closed subspace TC\\ of the total state space Ti is 



pJt) = Pp(t)P, 



(70) 



If the subsystem described by p\\(t) is an open system, i.e., if transitions from 
subspace TC\\ into TC Tt\\ (and vice versa) occur, then P cannot commute 
with the total Hamiltonian H. 

Now, in order to describe an n state system of the considered type, p 
of the form ( |55|) and a projector defining the subspace of the form (|27|), or 
another unitary one equivalent to it, should be chosen. It is easy to verify 
that for this P we have 

Po = PpoP, (71) 
so, in this case (see (1701) and (|69|)) 



P \\(t) = P P (t)P EE PU(t)P P0 PU + (t)P 



(72) 
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Using the identity (|9|) we have 




(73) 



{ dU\\(t) 
dt 



)p\\(t)+p\\(t)(i 



aunt) 



) 



(74) 



dt 



or, equivalently 



i-p^t) = H\\{t) p\\{t)-P\\{t) Htit) 



(75) 



(where H\\{t) is given by the identity ([13])), which is analogous to (|67D. 

6 Conclusions 

This paper deals with the operator form of the effective potential governing 
the time evolution in n-dimensional subspace of states. The general expres- 
sion for such an effective potential has been found in Section 2. Sections 3. 
and 4. dealt with the explicit form of such an operator for 2 and 3 dimen- 
sional cases. In Section 5. an application of the formalism developed in the 
previous sections to the density matrix has been suggested. 

The approach presented in this paper can be considered a natural exten- 
tion of the Wigner-Weisskopf approach to the single line width to more level 
subsystems which interact with the rest of the physical system. It has been 
shown that in the case of n level systems the WW approach may only be 
suitable if the PHP is n - fold degenerate, which of course is not always the 
case. 

The physical problem which is currently investigated with the use of si- 
milar methods is the neutral kaon complex and the possible violation of 
the CPT symmetry. This problem is obviously a 2 dimensional problem 
and can be researched with the use of the formalism developed in Section 
3. The standard approach to the problem developed in || uses the WW 
approximation to describe the time evolution of the K , K complex and 
proves to be quite a successful approximation of the physical reality. As noted 
at the end of Section 3., one of the conclusions which can be drawn here is 
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that hn = hii- This can be measured, and the parameter 5cpt ~ hu — ^22 
is used in tests of CPT conservation. However, if we want to retain the 
geometry of the problem (i.e. we do not want to reduce the problem to a 
one dimensional problem by assuming PHP degeneration) we will find that 
$cpt 7^ even under CPT conserved. For a more extensive discussion of 
this problem see ||, |j| [7|. 

The three dimensional case has not as yet been applied to describe an 
actual physical system and the possibilities of doing so will be investigated in 
future papers. One possiblility is to use the density matrix approach which 
has been proposed in Section 5., to the description of multi-level atomic 
transitions. Experiments designed to demonstrate the Quantum Zeno effect 
provide an example of such multi-level systems. For example Cook suggested 
an experiment which should demonstrate this effect on an induced transition 
in a single, trapped ion This experiment assumes the ion to have a 3 



level structure, and to describe it the density matrix approach is usually used 
(see for example HU). This gives us a possiblity to use the results obtained 
in Section 4.1 (Ai 7^ A2 7^ A 3 ) and Section 5. to construct a suitable equation 
for the reduced three dimensional density matrix. This, however, is beyond 
the scope of this paper and, as noted earlier, will be researched in future 
papers. 
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